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The smallness of the variation rate of the hamiltonian matrix elements compared to the (square 
of the) energy spectrum gap is usually believed to be the key parameter for a quantum adiabatic 
evolution. However it is only perturbatively valid for scaled timed hamiltonian and resonance 
processes as well as off resonance possible constructive Stiickelberg interference effects violate this 
usual condition for general hamiltionian. More general adiabatic condition and exact bounds for 
adiabatic quantum evolution are derived and studied in the framework of a two-level system. The 
usual criterion is restored for real two level hamiltonian with small number of monotonicity changes 
of the hamiltonian matrix elements and its derivative. 
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Adiabaticity is at the border between dynamics and 
statics. It has been introduced by Boltzmann in classi- 
cal mechanics and by Born and Fock in 1928 in Quantum 
Mechanics Q, Q , extended to the infinite dimensional set- 
ting by Kato (1950), studied as a geometrical holonomy 
evolution by Berry (1984), finally extended to degener- 
ate cases (without gapcondition) and to open quantum 
system more recently [3, 0] ■ The quantum adiabatic the- 
orem is usually used to derive approximate solutions of 
the Schrodinger equation and is strongly related to the 
(semi-)classical limit h — > of quantum mechanics [j| 
and to the Minimal work principle [f| for the Hamilto- 
nian H(t). The principle is simple: if a quantum system 
is prepared in an eigenstate |n(0)) of a "slowly" varying 
Hamiltonian it remains (without taking into account of 
the phase evolution) close to the instantaneous eigenstate 
\n(t)) of this Hamiltonian as time t goes on. The appli- 
cations range from two-level systems (nuclear magnetic 
resonance, atomic laser transitions, Born-Oppenheimer 
molecular adiabatic coupling, collisional processes ...) to 
quantum algorithms Q. 

"Usual" adiabatic conditions are (for all i £ [0,T]): 
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where the dot designs the time derivative and |m) are the 
instantaneous eigenstates for the energy eigenvalue E m (t) 
with E mn = fiWrnn = Em — E„ [271 . Some confusion oc- 
curs recently f^. floL fllL Il3i| because, this condition 
seems written for a general hamiltonian H(t). However, 
it has been studied by many different techniques (see for 
instance [Til fl5T|) but only for special types of hamilto- 
nian such as time scaling one H(t) = H(t/T) 2S\- Fur- 
thermore, even for such a time scaled hamiltonian, con- 
dition JIJ is not sufficient because it is only the leading 
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order term |l(il ITtI ] , in a time evolution T perturbation 
point of view, and more accurate conditions are needed 
to prove adiabatic evolution [l5j . 

The goal of this article is to derive general quantum 
adiabatic conditions for general hamiltonian. We start 
our study on a two level system example in order to study 
some possible violation of the usual adiabatic conditions. 
Afterwords, considering a more general type of N levels 
hamiltonians, we derive a general criterion for adiabatic- 
ity. Finally, the study of the interference during multiple 
passages allows us to precise the validity of the usual 
adiabatic condition. 

A quite general 2x2 hamiltonian matrix, written in 
the Pauli Matrix (a) basis, leads to the a spin 1/2 form 
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where luq = is the Larmor frequency, B is a rotating 
magnetic field with a polar angle 9, an azimuthal rotat- 
ing angular frequency tp = —ljl- Where the second form 
of the hamiltonian represents, in the rotating wave ap- 
proximation (RWA), a two level system coupled to an 
external (laser with angular frequency u>l for instance) 
field which is frequency detuned by Sq = c °s — ojl 
from the resonance and with a real Rabi frequency Qq = 
loo sin 8. For future developments we also define Qr, = 
sin 0-i0= |fi L |e iarg ^,(5 L = w L cos0-w o + ^ argfi L 
and Qji 



8 L . The eigenvectors e ^It); cor- 



responding to the eigenvalues =p^kj /2, are given by the 
columns of Re = I ,» c ° s . 2e „ e _ , v 2 sm 2 . e , + ]■ The 



evolution of the amplitudes and b + of the e |— ) 
and e l6+ 1 +) states are driven by the hamiltonian H = 
rIhRb - ihRliie: 
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9± is the "first order" choice 0±* = T5 Jq wo - ojl cos 9 
annulling the whole diagonal terms. 

Let us treat the (Schwinger 1937) example, where all 
the parameters loq, 6, <p = —lul are real and time indepen- 
dent. The evolution operator in the adiabatic | =p) basis 
U(t,0) = Rl(t)U(t,0)Re(0) (where U(t,0) is the evolu- 
tion operator in the diabatic basis) verifies ihU = HU 
and. with 9± — 0± \ is given by the matrix: 
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The adiabaticity (negligible off-diagonal terms in U) 
evolution is given by the following condition A' 2 ' = 
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yfio + ^o is the generalized Rabi frequency. Using = 



b\ this adiabatic condition can be written 
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which has to be compared with the "usual" adiabatic 
condition given by Eq. Q): 
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A(°> , A^ 1 ' , AW notations will be generally defined latter. 

Looking at the wi « wo and 9 very small resonant 
case (8l ~ So ~ 0), we see, in a simpler way than in Ref. 
0, and contrary to what is sometimes claimed [T2I Il3| , 
that Eq. is verified but not Eq. J3J. This fundamen- 
tal conclusion, based on a hamiltonian H(t) 7^ H(t/T) 
is still valid for the time scaling case H(t/T). Indeed, 
the Schwinger hamiltonian can be of the H(t/T) type 
if u)lT is taken to be constant, for instance by looking 
at the evolution after one period T = Tj, = 2-k/ljl de- 
pending on the u>l parameter value. Indeed, 2A^ = 



1 



sin 



+ 0(T" 2 ) 



1 



2A { 



+ 0(T~ 2 ) 



indi- 



cates, for instance if 9 is very small, why an evolution 
time Tl much longer than expected by the usual con- 
dition (T L > 
abatic evolution 



sin 9 



) can be needed to provide adi- 

The "usual" adiabatic conditions are 
then clearly not sufficient to provide adiabatic evolution 
even for H(t/T) hamiltonian type. 

To be more general let us now study a discrete, but 
possibly degenerate, hamiltonian with the state evolu- 
tion |tf(t)) = ELiM*)e <Mt) K*)> (N > 2). The 
phase 6 m = j m + a m is real but not necessary equals to 
the first order choice 9$ — Jq i(m\rh) — J* E m /h: geo- 
metrical phase (which is the Berry Phase for cyclic evo- 
lution) plus dynamical phase neither contains the (Pan- 
charatnam) phase arg{m(0)|m(i)}. To study the adi- 
abatic evolution we shall assume that \^f(t = 0)) = 



|n(0)) (i.e. b n (0) = 1). The evolution is adiabatic if 
1 - \(n(T)\m(T))\ = 1 - \b n {T )\ < 1 or equivalently if 
|||*)(*|-|n)(n||| = yT _ M«l[l3. 

The Schrodinger's equation leads for each m state to: 
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where 9 km = 9 k - 9 m . Using 9 m = 
J J and the norm inequality \/N 
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V^^/E™*, l&ml 2 > Em^lfcml we find the first 
(very restrictive) valid adiabatic condition for the inter- 
action time T: 



l-\K{T)\ < l-co8(^/N=lSl n T) < {N-lY-^T 2 (5) 

where Q n = max 4g [o,T] m ^ n l( n (0 1^(0) I- This condition 
is optimal because it is reached (see U) by the Schwinger 
N = 2 level system for 5 L = (ft„ = |ft L | = It 
illustrates the quantum Zcno effect: during a time much 
smaller than ^ n the system evolution is frozen. 

In order to find more useful adiabatic conditions we 
integrate by part Eq. J3J using (for k ^ m) A km = 
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It is now straightforward, with in = n, to look back 
to the standard adiabatic theorem with the time scaling 
t = sT. The evolution equation for |*(s)) = |*(t(s))), 
is then i££|*(s)) = ff(s)|*(s)) and the T -> +00 limit 
is similar to ?i — > 0. With jkm = E mn /h, we have (for 

(!h /I. _ A W _ ( m \k)e-I«"W-<"> l' h » 
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and the 



stationary phase theorem (saddle-point or steepest de- 
scent method) annuls, for T — > +00, the integrals in Eq. 
© leading to valid quantum adiabatic condition: 
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A comparison with Eq. Q indicates, as also shown by 
the two level model where \A^_\ = that a better 
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understanding of the o \Jp) term is in fact needed to have 
useful condition [l5j . 

We could now go back to the general H(i) case. 
sjl - \b n {T)\ 2 verifies yj\ - \b n (T)\ 2 < ^N~^\b_ with 
b- = max 4g r T i max m ^ n |6 m (i)|. Using Eq. (jfJJ) and 

Om = dm choice, it be bounded by 

2A + J\A'\ + (N -2)AflT 



If all A«, or Al n , 
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1 - (N - 2){A + JlA'l) - ((N - 1) + (JV - 2) 2 )AttT ' 



The typewriter style, such as (N — 1)AOT, indicates 
terms that can be annulled by using a better phase 
for 9 m namely the "second order" one dm = Orn + 

/„ f E ¥m e , (^+ fl -< 1) )4m(fc|™). The three impor- 
tant parameters are: 
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Where, Ai? = min^m E km is the energy spectrum 
gap. Another (better for large T) bound for b + (T) = 
mm tg[o,T] |frn(i)| = |^n(ir)| is obtained using m = n in 
Eq. © and the norm inequality: 

1 - |&n(*r)l ^ (N- 1)AQT 

+VF~^Ty / l - b\{A + VF" 7 ! J\A'\ + (N- 2)AttT) 

and a point fix study leads to 

1 - 6+ < 2(N - 1)AOT (7) 



+2(N - 1)(A + VF^T J\A'\ + (N - 2)AVLT) 2 . 
Finally one (not optimized) adiabatic condition is 

A + VNj\A'\ + (V^ + N-2)AflT<$: -j= (8) 
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We define two useful reals A km : 
(i) _ 



km 
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for the 9 km = 9 k ] } n choice , and A km = for the 



h m = Ofl choice where <y£[ = 7^ + E 
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When the hamiltonian H(t) is real in the canonical basis, 
the eigenstates \m) and (m\k) are reals and (m|m) = 
so, \A k 1 i\ = \A k ° m \. 



are monotonies in [0, T] 
Jo \ A km\ = \A km (T) - A fem (0)| and the condition © 
becomes simpler: A {1) + ^NA^flT < 1/AT or A (2) + 
y/N -2AWQT < 1/AT, where AW indicates that it 
should be calculated using the A km choice. For N = 2 
smallncss and monotonicity of A+l_ = ^j^- is equivalent 
to smallness and no more than one monotonicity change 
= > 0. Thus, a final general, 



simple and useful adiabatic condition is (for monotonies 
A km ) 

A {1) + \/N - 2A (1) fiT < 1/N. (9) 

It is even possible to refine the condition by dividing 
the interval [0, T] in smaller intervals where all A km are 
monotonies. A pcrturbativc point of view, neglectin g th e 
A^VtT term, has been used to derive similar results [18| . 

The N — 2 case is illustrative because it is the only 
one where a time independent adiabatic condition exists: 
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where M — 1 is the number of monotonicity change of 
, 00] . This generalize the Schwinger conditions 
For real hamiltonian the condition is 
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and becomes the usual adiabatic condition if M is small, 
for instance if the matrix elements b~o, fio of H and H 
have small number of monotonicity changes. This explain 
why the real dressed state hamiltonian, H = R\HR — 
ihRlk a = -§ (q° "| ), obtained from H in the rotat- 
ing frame (with the simple phase choice 8 + = 9- = 0) 
or simply by u>l = 0, have been luckily combined with 
the usual adiabatic theorem to describe several adiabatic 
evolutions such as, the RAP (Rapid Adiabatic Passage), 
the SCRAP (frequency or Stark-Chirped RAP) or the 
STIRAP (STImulated Raman Adiabatic Passage). 

However when real oscillatory terms are present the 
usual adiabatic condition is no more sufficient to pro- 
vide adiabatic evolution. As example we use the cycling 
hamiltonian flflL Eof . H = Hq with 5o(t) = acos(u>t) 
and a,LO,Qo are (positives to simplify) constants. It 
is relevant in many areas in physics: magnetic reso- 
nance, atomic collision, laser-atom interactions without 
the RWA and even localization by exchanging the pa- 
rameters So and f^o (hamiltonian R y HoR} y with R y = 
e iira v /4y rp^g wgak-coupling and large amplitude case 
a» S7o, to is simple because the non-adiabatic transition 
probability p% (so called single-passage or one-way tran- 
sition) is given by one of the simplest of the several ex- 
isting approximate formulas (Landau-Zener-Stuckelberg, 
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Rosen-Zener-Demkov, Nikitin, Zhu-Nakamura models, ... 
0H3) namely the Landau-Zener one: p\ ~ e~ 27r 3^ = 
e _ OT /(4AW(oo)) jj^J The M = 2 double-passage transition 
probability p2, which depends of a relative (Stiickelbcrg) 
phase of the wavefunction, p2 = 4pisin 2 (0) can be 
4 times higher than pi and the M (even) multiple pas- 
sage probability pu ~ Pi s1 c " s *q 9 can be M 2 times higher 
than pi . Here small lu value leads to the adiabatic limit 
Pi — > and with 8 ~ — ~ ir/2 we could have Pm ~ 1 
|22j . Interestingly enough, the reverse case, namely the 
diabatic limit (pi — > 1) can leads (for instance when a/w 
annul the Bcsscl Jo function) to the reverse phenomenum 
of adiabaticity created after multiple passages (j>m ~ 0) 
known as suppression of the tunneling, coherent destruc- 
tion of tunneling, dynamical localization or population 
trapping depending on the context 0, ^ • 

This two level example illustrate why monotonicity is 
require to avoid constructive interferences transforming 
an adiabatic (resp. diabatic) single passage in a fully 
diabatic (resp. adiabatic) transition after multiple pas- 
sages. The two level system with several crossings is 
very similar to the case of single crossing but with sev- 
eral levels leading to sum of dephased Landau-Dykhnc- 
Davis-Pechukas formulas [j^l l24j . Moreover, the tran- 
sition probability in a multilevel system is the product 
of several Landau-Dykhne type terms, corresponding to 
several successive transitions between pairs of levels |25f . 
However, several consecutive constructive interferences 
are exceptional and the generic most common case con- 
cern a system "complex enough" with small total proba- 
bility when the single crossing probability is small |26j . 

In conclusion, we have derived exact bounds for the 
evolution Eqs. JSJ, as well as general adiabaticity cri- 
terion Eqs. ©, (| 1 Of) - The key parameters for adiabatic- 
ity are the smallness and the small number of monotonic- 
ity change of A' 1 ' ~ — nr4S as well as a short evolution 



time (T » (N - 2) 



1 \\H\\ 
7 (i) AE 

3/2.1 ||H|| 



■). For real hamitonian 
is the spec- 
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the adiabatic (Pancharatnam) phase type 7 1 - 1 *' 
trum frequency gap and the usual adiabatic condition are 
restored if the matrix elements of H and H have small 
number of monotonicity changes in the two level (N = 2) 
case. The results presented here, and demonstrated for 
the discrete, but possibly degenerate case, might be use- 
ful for adiabatic quantum evolution and adiabatic quan- 
tum computation studies. Extension to the infinite di- 
mensional or non hermitian cases are some of the next 
steps needed to derive more universal quantum adiabatic 
conditions. 
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